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Magnetic control of ultra-cold ®Li and ^^"^Yb(^P 2 ) atom mixtures with Feshbach 
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We theoretically evaluate the feasibility to form magnetically-tunable Feshbach molecules in col¬ 
lisions between fermionic ®Li atoms and bosonic metastable ^^'^Yb(^P 2 ) atoms. In contrast to the 
well-studied alkali-metal atom collisions, collisions with meta-stable atoms are highly anisotropic. 

Our first-principle coupled-channel calculation of these collisions reveals the existence of broad 
Feshbach resonances due to the combined effect of anisotropic-molecular and atomic-hyperfine in¬ 
teractions. In order to fit our predictions to the specific positions of experimentally-observed broad 
resonance structures [1] we optimized the shape of the short-range potentials by direct least-square 
fitting. This allowed us to identify the dominant resonance by its leading angular momentum quan¬ 
tum numbers and describe the role of collisional anisotropy in the creation and broadening of this 
and other resonances. 

PACS numbers: 


The formation of ultracold molecules with a single un¬ 
paired electron allowing for coupling of the electron spin 
with the rotational angular momentum of the molecule 
has received increasing interest [ng. For example, the 
authors in Ref. [5] proposed that such spin-rotational 
splitting of molecular rotational states can make 
the long-range electric dipole-dipole interaction between 
molecules spin-dependent. Then, these molecules con¬ 
fined in two-dimensional optical lattices create a class of 
ultracold molecules that can realize spin-lattice models 
with unique topological properties. 

Ultracold molecules can be formed from an alkali- 
metal and an alkaline-earth-like atom in the ground or 
even long-lived metastable states. Experimental efforts 
toward production of these ground-state molecules are in 
their initial stages [I115HS]- The success of these experi¬ 
ments signihcantly depends on the realization of a two- 
step process, when a mixed quantum gas of alkali-metal 
and alkaline-earth-like atoms is hrst optically associated 
into weakly-bound molecules, which are then transferred 
into the rovibrational ground state. 

Photoassociation for homonuclear alkaline-earth-like 
molecules via so-called optical Feshbach tuning was pio¬ 
neered by [Mil]. This photoassociative tuning becomes 
possible due to the existence of long-lived excited molec¬ 
ular states near the narrow intercombination lines of the 
alkaline-earth-like atoms. Recent successful experiments 
[SHm showed that a two-photon optical Feshbach res¬ 
onance can be used to couple two colliding atoms to a 
vibrational level of the molecular ground state. The sup¬ 
pressed excited-state spontaneous decay makes efficient 
coherent molecular formation possible. 

An alternative method to form ultracold molecules 
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from ultracold atoms is by magneto-association. Mag¬ 
netic Feshbach resonances have had an enormous im¬ 
pact on the field of laser-cooled ultra-cold atoms and 
molecules [IIHII]. These resonances have allowed a tun¬ 
able, variable interaction strength by simply varying the 
strength of an external magnetic field (typically on the 
order of 100 G or 10 mT.) They are now widely used to 
investigate Efimov physics [22], create ultracold molec¬ 
ular gasses ona, and to simulate unique many-body 
phases [21]. So far, these resonances have been restricted 
to systems with ground-state alkali-metal and ground- 
state rare-earth species [201 . Only recently have 

experiments searching for resonances in collisions with 
electronic excited atoms sprung up. References [T] [30] 
detected the first resonances between ultra-cold ground 
state Yb or Li and excited Yb* in the metastable ^P 2 
state. 

The nature of Feshbach resonances in alkali and 
alkaline-earth-like atomic collisions significantly differs 
from that in alkali-metal atom collisions. In alkali met¬ 
als the hyperhne interaction between electron and nuclear 
spins gives sufficient complexity that leads to the appear¬ 
ance of Feshbach resonances. In alkaline-earth-like atoms 
such hyperhne interaction does not exist as the nuclear 
spin is decoupled from the electronic degree of freedom. 
Nevertheless, Refs. [SI [32] showed that magnetically tun¬ 
able Feshbach resonances in such systems can occur and 
are due to a weak i?-dependent behavior of the hyperhne 
coupling constant of the alkali-metal atom, where R is the 
interatomic separation between the atoms. These reso¬ 
nances are predicted to be narrow, on the order of mG, 
and most likely appear at large magnetic held strengths. 
Hence, they are difficult to observe and control. 

An interesting approach to induce and modify the po¬ 
sition and width of narrow Feshbach resonances was pro¬ 
posed in Ref. [53]. The authors show that intense non¬ 
resonant radiation can increase the resonance widths by 
three orders of magnitude up to a few Gauss. These ac- 
held-induced resonances can be used for the production 
of ultra-cold ground state molecules from the colliding 
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atoms. In addition, Refs. [MH55] suggested to apply an 
external dc-electric field with or without a magnetic field 
to control the interaction between ultra-cold atoms. 

Another promising way to observe broader and 
stronger magnetic Feshbach resonances and subject of 
this study is to consider interactions between one ground 
and one long-lived metastable atom. These metastable 
Feshbach resonances and its associated weakly-bound 
metastable molecule might be used to efficiently trans¬ 
fer colliding atoms to a vibrational level of the absolute 
molecular ground state. 

We report on widening the search for Feshbach res¬ 
onances in these exotic collisions by investigating het- 
eronuclear collisions of ground-state alkali-metal ®Li 
and long-lived metastable rare-earth ^^^Yb atoms. We 
demonstrate that the existence of magnetically-tunable 
Feshbach molecules for this system leads to formation 
of broad and strong resonances. Furthermore, we show 
that the origin of these Feshbach resonances is the result 
of a combination of exchange interactions and scattering 
anisotropies, where the interactions depend on the ori¬ 
entation of the interatomic axis relative to the magnetic 
field direction. The broad magnetic Feshbach molecules 
have non-zero orbital angular momentum 1. 

Our results confirm earlier predictions by Ref. m 
that resonances are strongly suppressed due to inelas¬ 
tic processes to energetically lower-lying ^Pi and ^Pq 
states. Nevertheless, a recent experiment [T] on the 
®Li-|-^^^Yb(^P 2 ) system has shown that resonances ex¬ 
ists. Here, we have performed least-square simulations 
on the shape of the molecular potentials and reproduced 
their resonance spectrum. In addition, we identified the 
strongest and dominant resonance in terms of its mixture 
of partial waves £. 

I. INTERACTION POTENTIALS 

For the determination of collisional properties of 
Li(^Si/ 2 )-l-Yb(^P 2 ) system we have (i) determined the 
short-range and long-range electronic interaction poten¬ 
tials V{R,t) and {ii) setup a close-coupling calculation 
that treats the hyperfine and Zeeman interaction, molec¬ 
ular rotation, magnetic dipole-dipole interaction, and the 
electronic interactions on equal footing. Here, R de¬ 
scribes the interatomic separation and orientation and 
T labels other quantum numbers that uniquely label 
the electronic potentials. We have calculated the non- 
relativistic doublet and ^H and quartet and 
“^H potentials, shown in Fig. using the configuration- 
interaction method within the MolPro package [55] . 

We used the aug-cc-pCVTZ basis set for Li [39] and 
chose a basis set constructed from the (I5s I4p I2d Ilf 
8g)/[8s 8p 7d 7f 5g] wave functions of M. Dolg and X. Cao 
[40l|4T] for Yb. The ytterbium basis relies on a relativis¬ 
tic pseudopotential that describes the inner orbitals up 
to the 3d^° shell. Only, the 2s valence electrons of Li 
and 4f^° and 6s^ valence electrons of Yb are correlated in 



R (units of Og) 

FIG. 1: (Color online) The ground and lowest excited non- 
relativistic potentials of the LiYb molecule as a function of 
internuclear separation R. Here ag is the Bohr radius of 
0.0529 nm. Solid blue and red curves are and ^H po¬ 
tentials, while dashed blue and red curves are and ^H 
potentials, respectively. Four metastable potentials dissociate 
to the Li(^Si/ 2 )+Yb(^P) limit and are used in our coupled- 
channel modeling. The green dashed horizontal line indicates 
the energy of a Feshbach molecule near the ^S-|-®P limit. 


the ab initio calculation. Our potentials agree well with 
the calculations of Ref. [52] • However, the potentials are 
not spectroscopically accurate and we have modified their 
shape to fit to the experimental loss rate spectrum. For 
R > 27ao, beyond the Le Roy radius where the elec¬ 
tron clouds of the atoms have negligible overlap, these 
electronic potentials are smoothly connected to the long- 
range isotropic and anisotropic van-der-Waals potentials 
determined from atomic polarizabilities [5]. In addition, 
we included the magnetic dipole-dipole interaction. 


II. RELATIVE STRENGTH OF VARIOUS 
INTERACTIONS 

Next, we find it instructive to analyze the relative 
strength of the long-range atom-atom interactions with 
the Zeeman, hyperfine, and relative rotational interac¬ 
tions. Figure]^ shows the two major anisotropic interac¬ 
tions, the magnetic dipole-dipole C^/R? (MDD) and the 
anisotropic dispersion ACq/R^ (AD) interaction, which 
can lead to the reorientation of the Li and Yb* angular 
momenta during a collision, as a function of R. These 
potentials as well as the isotropic dispersion C^jR^ (ID) 
interaction were drawn, assuming weighted-average val¬ 
ues of Cs and Cq, and ACe, based on our previous cal¬ 
culations [9]. In fact, Cg = (Cgs + 2C'6n)/3, where Cgs 
and CgE are the dispersion coefficient of the E and H 
potentials, respectively. The ACq coefficient is found 
from tensor analyses of the parallel and perpendicular 
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FIG. 2: Level splitting due to the dominant interactions 
between ®Li and Yb(^P) as a function of interatomic sepa¬ 
ration. Here Az is the Zeeman splitting at two magnetic 
fields, 5hs is hyperfine splitting of ground state ®Li, MDD 
is the magnetic dipole-dipole interaction Cs/R^, ID is the 
isotropic dispersion interaction CefR^, and AD is an es¬ 
timate of anisotropic dispersion interaction ACe/R^. Fi¬ 
nally, ROT is the rotational energy h?/{2^rR^)- We used 
Ce = 2694£;hao, ACe = 2425^0^, and C 3 = 6.6 x IQ-^Ehal 
Here Eh = 4.3610”^® J is the Hartree. 


components of the dynamic polarizability, as described 
in Section 2 of Ref. [43]. In addition, the molecule can 
rotate, here estimated by h'^/{2fXrR^) with reduced mass 
fir- Finally, hyperfine splitting (hs) of the ®Li atom, 
and the Zeeman interaction at magnetic field strengths 
B = 100 G and 400 G are shown. At much shorter sep¬ 
aration and larger energies, not shown in the figure, the 
exchange interaction due to the exponentially-increasing 
electron cloud overlap lifts the degeneracy of the double 
and quartet potentials and mixes the spin-orbit states of 
the Yb(^Pj) levels. 

The figure is most easily interpreted in the coordinate 
system and angular-momentum basis set with projection 
quantum numbers defined along the external magnetic 
field direction. In this coordinate system the rotational, 
spin-orbit, and hyperfine, and Zeeman interactions as 
well as the isotropic dispersion potential shift molecular 
levels, whereas the exchange interaction, the magnetic 
dipole-dipole interaction, and the anisotropic component 
of the dispersion potential lead to coupling between ro¬ 
tational, hyperfine, and Zeeman components. For dif¬ 
ferent interatomic separations different forces dominate. 
For example, when the curves for the magnetic dipole 
or anisotropic dispersion interaction cross the hyperfine, 
Zeeman and/or rotational energies spin flips combined 
with changes in the rotational state can occur. In fact, 
the anisotropic dispersion curve crosses the B =100 G 
and 400 G Zeeman curves near R = 50ao and 40ao, 
respectively, whereas the magnetic dipole-dipole curve 
crosses (5hs and Az at much shorter R, where chemical 
bonding or exchange interaction will also play an impor¬ 
tant role. 


III. CLOSE-COUPLING CALCULATION 

We have setup a close-coupling model for the scattering 
of fermionic ®Li and the metastable bosonic ^^^Yb iso¬ 
tope. The lithium atom is uniquely specified by electron 
spin SLi = 1/2 and nuclear spin Ili = 1. The metastable 
Yb* is specified by electron spin syb = 1 and electron 
orbital angular momentum lyb = 1- It is also convenient 
to define total atomic angular momenta fu = sli + tu 
and jVb = sVb + Wh for Li and Yb*, respectively. Their 
projections along the R-field direction are mLi and ruyb. 
The Hamiltonian for the relative motion of the two such 
atoms in a magnetic field B along the 0 direction is 

" = + ^ + 4) 

where the rotational Hamiltonian is /{2nR^), the 
spin-orbit interaction of Yb* is Hso = oso^Yb • syb, the 
hyperfine interaction of ®Li is H^f = OhfSLi • ?Li, and the 
Zeeman interaction Hz = (ffe.LiSLi,^ +ffAr,LHLi.z + ^Yb.z + 
ge,YhSYb,z)^J■BB^ with Bohr magneton ^lb and projection 
operators SLi,^ etc. The strengths aso and Ohf are taken 
from im US). The ge,Li, ge.Yb, and gN,Li are the elec¬ 
tronic and nuclear g-factors of Li and Yb* from [45| [46| . 
Finally, U{R) describes the magnetic dipole-dipole in¬ 
teraction and the non-relativistic electronic potentials, 
which are reexpressed in terms of (tensor) coupling op¬ 
erators between electron spins SLi and ^yb, and the an¬ 
gular momenta Zyb and i. For i? —^ 00 the interaction 
U{R) 0. 

There are four contributions to U{R): a spin- 
independent isotropic potential Viso{R), which is pro¬ 
portional to an attractive 1/i?® potential for large sep¬ 
arations, a short-range isotropic exchange interaction 
G sV&]oo = -V^’J^{R)sLi ■ SYb/VS, which 
splits doublet from quartet E and H potentials and falls of 
exponentially for large R, and an anisotropic quadrupole- 
like interaction t 4 ni(R)[G 2 (.R) G [lyb G ^y 6 ] 2 ]oo, which 
lifts the degeneracy of the E and H potentials and 
I4ni(R) oc 1/i?® for large R. Here, Ckq{R) is a spher¬ 
ical harmonic. The four interaction strengths Viso{R), 
14 xc^(i?), and Vani(R) are constructed such that in the 
body-fixed frame with projections along the internuclear 
axis the sum of the interactions closely reproduces our 
four non-relativistic potentials. 

The Hamiltonian is constructed in the atomic basis 
or channels 4>)\aLi)\o:Yb), where Yirai{d,4') is a 

spherical harmonic and angles 9 and (j) give the orienta¬ 
tion of the internuclear axis relative to the magnetic field 
direction. Hence, the rotational Hamiltonian is diago¬ 
nal in this basis. The kets |aLi) are eigen states of the 
atomic Zeeman plus hyperfine Hamiltonian of ®Li. For 
fields B > 100 G, of interest in this paper and where 
the Paschen-Back limit holds, we denote states by OLi = 
ms,mi, where iris and rrii are the projections of sli and 
?Li, respectively (and mu = ms+mi). The kets |ayb) are 
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FIG. 3: Adiabatic potentials of colliding ®Li(^Si/ 2 ) and 
^^^Yb(^P) near the ^P 2 dissociation limit at a magnetic field 
of B = 300 G and Mtot = —1/2 as a function of interatomic 
separation R. Channels with even partial wave £ up to 6 are 
included. Curves approaching the same dissociation limit cor¬ 
respond to different partial waves. The collision starts in the 
s-wave mLi = 1/2 and jVb = 2, rriYb = — 1 channel indicated 
by the dashed horizontal line. A resonance occnr when the 
energy a bonnd state of a closed channel, which are schemat¬ 
ically shown as dark-green horizontal solid lines close to the 
entrance-channel energy. Here k = 1.38065 • 10“^® J/K is the 
Boltzmann constant. 


eigen states of the atomic spin-orbit and Zeeman Hamil¬ 
tonian of in the term. As the spin-orbit in¬ 

teraction is orders of magnitude larger than the Zeeman 
interaction we denote the states by ayb = jYbWYb- Cou¬ 
pling between the basis states is due to U{R,t). The 
Hamiltonian conserves Mtot = WLi + WYb + and only 
even (odd) ^ are coupled. 

In the experiment of Ref. [1] the ®Li atoms are pre¬ 
pared in the energetically-lowest hyperhne state with 
mui = 1/2 while the Yb atoms are prepared in the ^P 2 
(bvbWYb) = |2, — 1)) sublevel. For ultra-cold gasses with 
temperatures near 1 /rK we can focus on so-called s- 
wave or £ = 0 scattering. Theoretically, this corresponds 
to calculations including channels with Mtot = —1/2 
and even partial waves £. We will also present results 
for ultra-cold scattering with Yb* prepared in the state 
IjYbWYb) = |2, —2). This corresponds to calculations for 
Mtot = -3/2. 

Figure shows a sample of the medium- to long-range 
adiabatic potentials, obtained by diagonalizing H exclud¬ 
ing the radial kinetic energy operator, as a function of 
interatomic separation R near the Yb*(^P 2 ) dissociation 
limit at a magnetic field of B = 300 G. Only poten¬ 
tials dissociating to the fu = 1/2 or 3/2 and jyb = 2, 
mYb = 0, —1,—2 limits are visible in the figure. Inelas¬ 
tic collisional losses to energetically lower-lying Yb*(^P, 
jYb = 0,1) states do nevertheless exist. Feshbach reso¬ 
nances occur in channels with a dissociation energy above 
that of the entrance channel. 

Figure shows the experimental inelastic rate coef- 



FIG. 4: (Color online) Upper panel: Loss rate coefficient in 
the ®Li(ms = — 1/2, rrii = l)+^^‘^Yb(®P 2 ,/ = 2,m = —1) col¬ 
lision as a function of magnetic field B. The experimental 
data (blue dots with one standard-deviation uncertainties) is 
obtained at a temperature of T = 1.6 ^K. The theoretical sim¬ 
ulation (solid red curve) with optimized short-range potentials 
is obtained for the same temperature. Lower panel: Scatter¬ 
ing length as a function of B at collision energy E/k = 1.6 
pK. The short-range potentials are as the same as for the rate 
in upper panel. Two dashed line indicates the van der Waals 
length with average value of G6=2694 a.u.. 


ficient [T] at a temperature of 1.6 in the collision 
between ®Li and atoms and our best-fit theo¬ 

retical inelastic rate coefficient at a collision energy of 
E/k = 1.6^K. We include channels with Mtot = —1/2 
and even £ up to 8 leading to a close-coupling calculation 
with 99 channels. A thermal average of the theoretical 
inelastic rate coefficient is not required as we are in the 
Wigner threshold limit where this rate is independent of 
collision energy. We observe at least one clear resonance 
aX B = 450 G and possibly three weaker resonances. 
From calculations that include fewer partial waves we 
find that the resonances can not be labeled by a single 
partial wave. Their locations shift signihcantly and only 
converge to within a few Gauss when £ = 8 channels are 
included. Including channels with £ < 10 shifts resonance 
by no more than 0.1 G. This is well within the width of 
the resonances for both Figs. 4 and 5. Gonsequently, 
we conclude that it is sufficient to include channels up to 
£ < 8 . 

Our bound-state calculation allows us to identify the 
dominant closed channel that couples to the s-wave en¬ 
trance channel at H = 450 G. We identify this chan¬ 
nel with the £=4 partial wave and projection to^= 1. 
This coupling is due to anisotropic interactions, such 
as the anisotropic dispersion forces at large separations. 
Moreover, this channel dissociates to the atomic limit 
®Li(TOs=l/2, mi=-l) -\- ^’^"‘Yb(^P, j=2, m=-l), which lies 
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above the entrance channel (see Fig. 3). 

Figure]^ also provides evidence that the inelastic rate 
of the strong resonance is close to the universal rate, as 
shown by the dashed line in Fig. This is a clear indi¬ 
cation that the collisional system is close to the universal 
regime for losses on the dominant resonance and lies well 
below for other magnetic fields 01117]. 

The agreement between the experimental [Tj and our 
theoretical spectrum is only obtained by allowing the co¬ 
efficients Viso{R), and Vs,ni{R) to vary such that 

the depth of the four non-relativistic potential curves 
is changed. We can not exclude the existence of other 
shapes of potentials which will lead to a loss rates that 
is consistent with the experimental data. 

The lower panel of Fig. shows the scattering length 
a as a function of the magnetic field strength for 
the collisions between ®Li(ms = —l/2,TOi = 1) and 
^^^Yb(^P 2 ,jYb = 2,TOYb = —1) atoms at a collisional 
energy of 1.6 /iK. Here the scattering length is defined 
through 5'eiastic = exp(-2*fc[a -ib]), where S'eiastic is the 
elastic S-matrix element for the s-wave scattering chan¬ 
nel, wavenumber k is given hy E = and pos¬ 

itive length b can be related to the inelastic loss rate co¬ 
efficient. The resonance features of a correspond to those 
for the inelastic rate in Fig. The Feshbach resonance 
centered at 450 G has a well resolved Fano profile with a 
background value close to 75ao. 

We have also studied the role of the Yb spin-orbit 
interaction and anisotropy of the electronic potentials 
on the Feshbach resonance structure. Our analyses in¬ 
dicate that for most magnetic fields the resonance fea¬ 
tures are broadened by the losses energetically lower col¬ 
lision channels with non-zero orbital angular momen¬ 
tum i. This confirmes our predictions about the im¬ 
portance of anisotropic coupling between the short-range 
potentials. Moreover, we also observe that, except for 
small fields, more than one half of the loss goes into 
^P 2 ,j = 2,m = —2 channels. In addition, a significant 
fraction of the population ends up in the Yb(^Pi) spin- 
orbit state (see Fig.4(b) in Ref. [j). 

IV. CONCLUSION 

We discussed advances in the understanding of scatter¬ 
ing properties of high-spin open-shell atomic systems. In 
particular, our attention was directed towards collisions 
between fermionic lithium atoms in their absolute ground 
state with bosonic ytterbium in its meta-stable ^P 2 state. 
Both species, either with their bosonic or fermionic iso¬ 
topes, have been successfully cooled to quantum degener¬ 
acy and confined in optical dipole traps or optical lattices, 
allowing the study of their collisions at the quantum level. 
We observed magnetic Feshbach resonances in collisions 
with the meta-stable ^^"^Yb atoms even though this atom 


was prepared in the myh = — 1 magnetic sub level, which 
is not the energetically-lowest ^P 2 Zeeman state. With 
couple-channels calculations we have elucidated the role 



FIG. 5: (Color online) Predicted theoretical loss rate coef¬ 
ficient in the ®Li(ms = 1/2,= — l)-|-^^‘^Yb(®P 2 ,/Yb = 
2,mYb = —2) collision with Mtot = —3/2 and collision en¬ 
ergy of E/k = 1.6 fiK as a function of magnetic field B. The 
simnlation is performed with optimized short-range potentials 
as in Fig. 

of anisotropies in the creation of these resonances and by 
optimizing the shape of the short-range potentials repro¬ 
duced the measured loss rate for magnetic fields up to 
500 G. 

Initially, collisions between ground state Li and 
metastable Yb in the ^P 2 state seemed to be a good sys¬ 
tem in which to form strong Feshbach resonances, since 
the Yb life time is 15 s and interactions with Li atoms is 
highly anisotropic. Our calculations, however, show that 
when Yb atoms are not prepared in the energetically- 
lowest magnetic sublevel the resonances have large in¬ 
elastic losses (on the order of 10“^° cm^/s). This large 
rate coefficient might be an obstacle to use them in stim¬ 
ulated Raman adiabatic passage schemes to create the 
ground state molecules. On the other hand, we also found 
that inelastic loss for the resonances in collisions with the 
energetically-lowest magnetic sublevel of Yb is an order 
of magnitude smaller. These resonances if experimen¬ 
tally confirmed can be promising candidates for Raman 
transfer. In Fig. we present a prediction for the Fes¬ 
hbach spectrum for the collision between ®Li and ^^^Yb 
(^P 2 ,j = 2,m = —2) based on the optimized potentials. 
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